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Classical and quantum dynamics of two distinct BPS monopoles in the case of non- 
aligned Higgs fields are studied on the basis of the recently determined low energy 
effective theory. Despite the presence of a specific potential together with a kinetic 
term provided by the metric of a Taub-NUT manifold, an 0(4) or 0(3,1) symmetry of 
the system allows for a group theoretical derivation of the bound-state spectrum and 
the scattering cross section. 
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Recently, much attention has been given to the low energy interaction of BPS monopoles in the 
case that a semi-simple gauge group G of rank k breaks to its maximal torus U(l) k by one or more 
Higgs fields in the adjoint representation. In the simplest non-trivial case of the rank two and when 
the Higgs expectation values have only one independent component, the relative dynamics of two 
distinct fundamental BPS monopoles can be described by the kinetic Lagrangian provided by the 
metric of a Taub-NUT manifold [jj], (2|. Here, relevant variables consist of the relative position r 
and the relative phase ip G [0,47r], and the Lagrangian takes the form 

^kin = f | (l + "f) r 2 + rl (l + r -f) 1 $ + w(r) • r) 2 j , (1) 

where /i is the reduced mass, and w(r) is a Dirac monopole potential satisfying V x w = — r/r 3 . 
The positive parameter tq of length dimension provides the length scale of the problem. The metric 
of the Taub-NUT space is chosen so that Lki n = fig^i^x" /2 with the coordinates = (r,roip). 
On the other hand, if the Higgs expectation values have more independent components, it has been 
shown recently || that the low energy Lagrangian acquires an additional attractive potential term. 
(For more than two distinct monopoles, there can be more than one term.) That is, the Lagrangian 
is of the form 

L = L kin - U(r) (2) 

with 

where a > is a constant related to the degree of the vacuum misalignment. 

In this paper we will present the analysis on classical and quantum dynamics of the system defined 
by the Lagrangian @ and exhibit physical consequences due to the potential (||) . 

With U(r) = and the parameter ro taken to be negative, our system reduces to the previously 
considered case in Ref. [§, [|, || for an approximate description of the relative dynamics concerning 
two identical BPS monopoles. In this case the system was found to possess the 0(4)/0(3,l) 
dynamical symmetry analogous to the one in the well-known Coulomb problem, and this in turn 
allowed an algebraic derivation of the dyonic bound state spectrum and the two monopole scattering 
cross section as shown in Ref. [|], [|, 0. This analysis requires no essential modification for our case, 
except that no dyonic bound states become available when the potential U(r) is absent. In our 
case the relative charge leads to a repulsion instead of an attraction between two BPS monopoles. 
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What changes would the presence of the potential @ in the case of non-aligned Higgs fields bring? 
What is remarkable about our system is that the 0(4)/0(3,l) dynamical symmetry operates even 
in the presence of this potential and so we can determine the effect due to the potential exactly. 
Moreover, we will see that at the level of a time-independent Schrodinger equation, our system 
shares the same form as the Zwanziger system QJ-a model involving a suitably combined monopole+ 
Coulomb +l/r 2 potential with 0(4)/0(3,l) symmetry. Based on this observation, we can deduce 
the exact spectrum of dyonic bound states (which exist only in the non-aligned Higgs vacuum), 
and also the scattering cross section for two distinct BPS monopoles in the general case. 



We will study the classical dynamics first. Suppose that we have not yet posited the precise 
form of a radial potential U(r) in Eq. (|2|). The Lagrangian has two obvious conserved quantities 
associated with the cyclic variable ip and the time t, viz., 



E 



fd + - 

2 r 



Q = V r o( 1 + — ) 1 (^ + w(r) -r), 

r 

r r 
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+ U(r), 



(4) 
(5) 



corresponding to the relative electric charge and energy, respectively. The mechanical three mo- 
mentum 

,r = M(l + -)r (6) 

r 

satisfies the equation of motion 

„2 \ 



7T 



fJ>rp r / . 2 
2 r 3 F 



r x r 



H 2 rl 



r U'{r). 



(7) 



There is also the conserved angular momentum, 



J = r x 7r + qr, 



(8) 



and so the position vector r moves on a cone defined by 



J • r = q. 



0) 



At this point, let us further demand, as a restriction on the potential U(r), that all bound trajec- 
tories in the system be closed orbits. Such is the case if there is an additional conserved quantity 
of the Runge-Lenz type, K = 7r x J — (5r with a constant [3. A simple calculation using Eq. ([?]) 
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K = 7T x J - nr [ E - I r. (10) 



determines the potential to be of the form (0) up to a possible constant shift. The precise form of 
the conserved Runge-Lenz vector is also fixed as 

t* r o, 

(Essentially, this corresponds to a generalization of the Bertrand theorem || to the case with the 
Taub-NUT kinetic term.) Our potential thus assumes a very special status-it corresponds to the 
general attractive potential which admits both the angular momentum and the Runge-Lenz vector 
as conserved quantities. 

One may wonder whether there exists any additional modification of the Taub-NUT dynamics 
without jeopardizing the existence of a conserved Runge-Lenz type vector. One such modification [EJ 
is to change the Taub-NUT metric to 

ds 2 = ^L dY 2 + ar + br 2 + 2 

r 1 + cr + ar z 

For a particle moving in the space with this metric, there still exists a conserved Runge-Lenz 
vector |J. However, the metric is always singular somewhere near the origin. This modification 
is rather ad hoc and does not have any good physics motivation. This metric also has nonzero 
Riemann tensor R^ u and so is not hyperkahler. Further exploration of this modification was 



studied in Ref. j|, 10 1. Another possibility, which is somewhat novel, is to add a term linear in 
velocity, 



^linear = m + w(r) ■ r). (12) 

r 

This term can be absorbed by shifting ip to ip + ci/(/Wg) and a 2 — > o 2 + c 2 in the Lagrangian 
(0). Thus the physics does not change much due to this modification. Thus, we will focus on our 
Lagrangian (||). 



Even in the presence of the potential (||), a complete description of the motion can now be given 
on the basis of the conservation laws. From two conserved quantities J and K, we get the conserved 
quantities 

J • K = -nr q{E - -^), (13) 
K 2 = 2^(J 2 - <f)(E - 4±#) + Ao 2 (£ " — 2 f- (14) 
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Defining a conserved vector 

N = qK + fir ( E J J , (15) 

we also find that 

Nt = q{3 2 -q 2 ). (16) 

From this equation, it follows that the motion is confined on a plane perpendicular to the constant 
vector N. Thus the trajectories belong to a plane intersecting the cone, i.e., correspond to conic 
sections. The detailed orbit characteristics can be inferred from the energy equation (||). From 
Eq. (||), the effective potential for a given q is 

^ = 4 (1 + ? )+ 4tt?' < 17 > 

For \q\ > a, U e fi(r) is a monotonically decreasing function of r and so there is no bound state. 
In this case, the energy is bounded by the minimum of the effective potential U(oo), i.e., E > 
(q 2 + a 2 )/(2/xrg). When the equality holds for this bound, the trajectories are parabolic orbits, 
and otherwise we find the hyperbolic orbits. On the other hand, if the relative electric charge 
satisfies \q\ < a, the potential U e d(r) assumes the minimum value E^, ps = a kl/(/ ir o) a * finite radius 
r bps = r o/( a /\l\ ~ !)• Thus, with \q\ < o, we have the following possibilities, 



E > 7r^z{q 2 + a 2 ): (unbounded) hyperbolic orbits, 

E = 2^5- (q 2 + a?): parabolic orbit, 

2^1 {q 2 + o?) > E > Eb ps : (bounded) elliptic orbits, 

E = Eb ps : a static configuration with r = 0, r = rb ps and J = qr. 



The last case where E = -Ebps can be regarded as the classical 1/4 BPS configuration |Tl| , H], 

For the scattering orbits, the two conserved quantities J and K can also be utilized in determining 
how the scattering angle O depends on the impact parameter b, the initial relative speed v (or the 
initial energy E), and the relative charge q. Since the detailed arguments on this can be found 
from Ref. Q, we shall here give only the formula pertaining to our case: 
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where we have used the relation 



1 



a 2 + q 2 \ 2 4a 2 q 2 



1 



(19) 



E = ^fiv 2 + 77^-2 + a 2 ) 

"*"o 



2 r ~~ ' 2 fir 2 



(20) 



for the initial variables. This leads to the Rutherford-type differential cross section 

da 
dtt 




2J2 



4a z q 
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16 | \ /(-/•,-<" / fi^r^v 
Observe that when o vanishes, this reduces to the old result |Q, ||, ||]. 



(21) 



(22) 



Let us now turn to quantum dynamics. For the classical Lagrangian (2), the canonical momentum 
conjugate to r is given by 

fir, 

(23) 



P = -^7- = 7T + <?w(r), 



Or 



and the Hamiltonian by 



H 



1 

2/i" 



>,, r Tn 



(24) 



2/i(l + r /r) , ^ , v 

In the quantum theory, dynamical variables become operators and so, for an expression involving 
operators like our Hamiltonian, an appropriate ordering must be specified. We will here follow the 
standard procedure developed for a general nonlinear system |l2|]. The variables x^ = (r, roip) and 
their conjugates p^ = (p, q/ro) are taken to be hermitian operators with respect to the Taub-NUT 
volume measure ^fgd 4 x = (1 + ro/r)d 4 x with g = det(g fll/ ), and satisfy the basic commutation 
relations = 0,\Pfi,Pu] = and [x^,p u ] = ihStf. (From now on, we put % = 1.) In this 

representation, the canonical momentum operator is 



Pn = -i(dn + -g l d tx g). 



The Hamiltonian operator is then identified with [12] 



H 



2/i 



g 1 Pny/ggTpvQ 



-1/4 



+ 



1 



^r 2 1 + =? 
a 2 



(25) 

(26) 
(27) 



The relative charge operator q = r^p^ = —id/dip commutes with the Hamiltonian and has eigen- 
values s with s = 0, ±1/2, ±1, .... If one also introduces a non-hermitian operator 



7T = — iSJ — qw(r) 



which satisfies the commutation relations 

[iri,x j ] = -iSj, [7v,ip] = iw, [n^iTj] = - i 



x k 



(28) 



(29) 



it is a matter of lengthy computations || to show that the Hamiltonian (27) can be recast precisely 
into the form (24), now with the operator ordering taken into consideration. The Hamiltonian 
operator in Eq. p3) is hermitian with respect to the Taub-NUT volume measure. 



We are interested in the solution of the time-independent Schrodinger equation, Hty = E^> . For 
the eigenstates with q = s, we use the Hamiltonian (^4[) and factor out the trivial phase e iq ^ from 
the wave function ^ to obtain 



2/j, r zfir 



77 - - + TTTZa ~ 8 r * = °> ( 30 ) 



where two new constants a and £ are given by 



-^2), S-E-- — 



a = r (E-^), S = E--^(a 2 + q 2 ). (31) 



Disregarding the above relations, one may look upon Eq. pQ ) as the Schrodinger equation of a 
flat-space system, TC^f = £^f, with the Hamiltonian, 

n ~Y^ "7 + 2^- (32) 

This Hamiltonian is hermitian with respect to the usual Euclidean three space volume measure 
and contains, aside from the magnetic monopole and Coulomb-type potentials, a 1/r 2 repulsive 
potential of a definite strength. The very system was first considered by Zwanziger M|; Eq. ( ]32| ) 
describes the relative dynamics in his exactly soluble nonrelativistic model of two dyons. It enjoys 
the 0(4:)/ 0(3, 1) symmetry. Our system also enjoys the 0(4)/ 0(3, 1) symmetry, which was noticed 
earlier without the potential @ 0, n] . In our case the symmetry generators are quantum operators 
corresponding to the angular momentum and the Runge-Lenz vector: 

J = r x 7T + qr , (33) 

1 Q 1 
K = -(tt x J - J x tt) - nr r(H - -M, (34) 
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where H is our full Hamiltonian operator (24). (Note that the order of various operators is im- 
portant now.) Then by a bit of lengthy algebra, one can verify that they satisfy the following 
relations: 

[J,H] = [K,H]=0, (35) 

[Ji,Jj] = ie ijk J k , (36) 

[Ji,Kj] = ie ijk K k , (37) 

2 2 

[K u Kj] = -ie ijk J k {2fxH - ^-4^). (38) 



They also satisfy 



Q 2 

K-J = 3-K = -^r q{H--^), (39) 



K 2 - frl(H - 4) 2 = 2MJ 2 " <? + h 2 )(H - q ^±^). (40) 



With the explicit % correction in the quantum expression term, one can compare these equations 
with their classical counterparts ( |i~3| ) and fll4|). For the Zwanziger system, the corresponding 
0(4)/0(3, 1) generators are made of J and 

K = -(7T x J — J x 7r) — fj,ar, (41) 

which is obtained from K by replacing the operator rg(H — q 2 / 'fir^) with the parameter a of 
Eq. (|3l|). They commute with 7i and satisfy the similar algebra with [Ki,Kj] = —ieij k J k (2fi7i), 
and so on. The analogue of Eqs. (|39|) and (^0| ) for these operators are 

K J = J K = -ptaq (42) 
K 2 - [i 2 a 2 = 2fi(3 2 -q 2 + h 2 )H (43) 



As in the Coulomb system [13], the bound state spectrum and the scattering cross section for the 
Zwanziger system can be found exactly with the help of an appropriate group representation theory. 
Then, they may be recast using the connection (|3l|), to obtain the corresponding results for our 
original system with the Hamiltonian (24). Of course, it is possible to use the group representation 



theory directly with our system also. But the procedure just mentioned is simpler. 
For a > 0, the Zwanziger system allows bound states with negative energy £ < 0. In this case, 



the operator M = K/y2ju|£| and J satisfy an 0(4) algebra. Introduce A = (J + M)/2 and 



B = ( J — M)/2 as the generators for two 0(3)'s of 0(4) = 0(3) x 0(3). Then an energy eigenstate 
may be taken as \a, 0,3, b, 63 > with integers or half-integers for 0,03 and 6,63, satisfying 

A 2 ^ = a(a + 1)^, A 3 ^ = a 3 $), (44) 
B 2 * = b(b + 1)*, 5 3 tf = 63*. (45) 

Then, making a judicious use of the identities ( |42|) and ([43]), it is not difficult jjj to derive the 
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exact bound state spectrum £ = — , where n is the principal quantum number restricted to take 
values n = |s| + l,|s|+2,... if the eigenvalue s of the operator q assumes half-integers or integers. 
This bound state energy levels are highly degenerate, since, for given s and n, the quantum number 
j associated with J 2 = (A + B) 2 can run from |s| to n — 1. We can now translate these findings 
into those for our original system. First of all, the conditions to allow bound states, a > and 
£ < 0, translate into 

^ <E< 2^ ( " 2 + S2) ' (46) 
and so our system can have bound states only with \s\ < a as in the classical case. On the other 
hand, inserting the formula £ = —t^z to Eq. ([H]), we obtain a quadratic equation for E, 

-^(E-— 2 ) 2 = E- -^(a 2 + , 2 ). (47) 

Solving this equation, we immediately obtain the desired bound state spectrum for our system: 

E = \-{n 2 - s 2 ) + nVn 2 - s 2 + a 2 } (48) 

with |s| < a and n > \s\ + 1. The degeneracy of each energy eigenstates is n 2 — s 2 as j runs from 
\s\ to n — 1. 



The case with £ > corresponds to a scattering state. We may put £ = iiv 2 /2 with the initial 
speed v. The operators M and J now satisfy the 0(3,1) algebra. As in Ref. 0, the group theory 
can be utilized to derive the exact quantum mechanical scattering cross section formula for the 
Zwanziger system, 

Again we translate this result into our case by using the relation (^) and q = s and see that this 
coincides with the classical cross section formula (|22"|), making the classical result to be quantum 
mechanically exact. 
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For a given value of s, the bound state energy spectrum in Eq. ( (48|) is an increasing function of 
n, satisfying the inequality 

E(n = \s\ + 1) < E < E(n = oo) = — Ur(a 2 + s 2 ). (50) 

Especially, the ground state energy E(n = \s\ + 1) with degeneracy 2|s| + 1 is greater than the 
classical BPS energy Eb ps = a|s|/(//rg), which is identical to E(n = \s\) if n could take the value \s\. 
The difference between the ground state energy and the classical minimum is due to the zero point 
energy of the bosonic theory. It would be interesting to formulate the most general super symmetric 



Lagrangian M, 14, rL5j for our system and find BPS and non-BPS bound states. There would be 



no such energy difference in supersymmetric vacua. 



Finally, there exists an even larger conformal symmetry group 0(4, 2)/0(5, 1) for the Taub-NUT 
system This can be easily generalized to our system, which we do not discuss here. They 

will again play a role in the supersymmetric extension of our work. 
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